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Quasinormal modes are the counterparts in open systems of normal modes in conservative sys- 
tems; defined by outgoing-wave boundary conditions, they have complex eigenvalues u. The con- 
ditions are studied for a system to have a supersymmetric (SUSY) partner with the same complex 
quasinormal-mode spectrum (or the same except for one eigenvalue). The discussion naturally 
includes total-transmission modes as well (incoming at one extreme and outgoing at the other). 
Several types of SUSY transformations emerge, and each is illustrated with examples, including the 
transformation among different Poschl-Teller potentials and the well-known identity in spectrum 
between the two parity sectors of linearized gravitational waves propagating on a Schwarzschild 
background. In contrast to the case of normal modes, there may be multiple essentially isospectral 
partners, each missing a different state. The SUSY transformation preserves orthonormality under 
a bilinear map which is the analog of the usual inner product for conservative systems. SUSY 
transformations can lead to doubled quasinormal and total-transmission modes; this phenomenon 
is analysed and illustrated. The existence or otherwise of SUSY partners is also relevant to the 
question of inversion: axe open wave systems uniquely determined by their complex spectra? 
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I. INTRODUCTION 

A. Supersymmetric quantum mechanics 

In this paper we consider supersymmetry (SUSY) in the one-dimensional Klein-Gordon equation (KGE) 

[d^,-dl + V{x)\<P{x,t) = Q , (1.1) 
and especially in the corresponding eigenvalue problem 

HM^) ^ Ljl(j)r,{x) , (1.2) 

where 

H = -dl + V{x) . (1.3) 
The boundary conditions will be specified later. In so far as the interest centers on the time-independent problem 



(1.2) and the spectrum, the case of the Schrodinger equation, to which reference is usually made, is included if the 
eigenvalues are simply relabeled by ui^ w. 
If there exists another system described by 

H^^dl + V{x) , (1.4) 

such that H and H have the same spectrum (or the same spectrum apart from one state) , and moreover if the states 
in the two systems are related by 

^ix)=A(j)ix) , (1.5) 

where 

A = d, + W{x) , 
-A^=d,-W{x) , (1.6) 
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then the two systems are said to be SUSY partners [Q. The term "super" originates from boson-fermion symmetry 
in field theory, but here simply refers to a relationship between the states of two different systems, whereas symmetry 
relates states of the sa me system. SUSY for quantum mechanics is usually discussed for normal modes (NMs), i.e., 
solutions to (1.1) and (1.2) which vanish at spatial infinity and which are square- integrable The results for NMs 
are well known ^ and will be covered in the discussion below. 
In order for (1.5) to preserve the spectrum, one needs 



AH = HA . 

Our task is simply to analyse this condition when applied to open systems. 



(1.7) 



B. Quasinormal modes and total-transmission modes 

In open systems, waves are not confined, but can escape to infinity: acoustic waves escape from a musical instru- 
ment; electromagnetic waves escape from a laser; and linearized gravitational waves propagating on a Schwarzschild 
background escape to infinity and into the horizon. These systems are often described (e.g., in the case of gravitational 
waves [^) by the KGE (LI), or (e.g., in the case of optics §) by the wave equation 



[p{x)d^-dl]<j,ix,t)^0 , (1.8) 

to which the KGE can be transformed . This transformation has certain subtleties when there are bound states in 
the KGE; see Appendix ^. 

We assume for the moment that V (or in the case of the wave equation p — 1) has finite support. This assumption 
is natural for describing a system in the finite parts of space, surrounded by a trivial medium such as vacuum, and 
it leads to a simple formalism. However, some of the superpartners of finitely supported potentials, as well as the 
examples of Sections ^ and VI, prompt the study of SUSY under the weaker condition V{x) ^ as a; ^ ±oo. The 



analysis then becomes considerably more subtle, and is therefore relegated to Appendix |_ 

We shall consider first of all states satisfying the so-called outgoing-wave condition, which is most easily stated for 
an eigenfunction of frequency u: 

(j){x) cx e=^'"^ , x ^ ±00 , (1.9) 

or equivalently 

—^^±iLO , X^±O0 . 1.10) 

cj){x) 

Under these boundary conditions, discrete eigenvalues fall into two classes. First, there could be bound states or 
NMs 0; because V vanishes at infinity, these must (from the Schrodinger point of view) have a negative energy 



and hence uj is purely imaginary. According to the condition (l.£), we put Imw > 0. Second, there could be 
quasinormal modes (QNMs) with complex eigenvalues uJ^ . Because these waves (in contrast to the bound states) are 
genuinely outgoing, probability (or energy) in the finite parts of space decreases, so Imw < We here assume that 
the loss is only due to the escape of the waves, i.e., only through the boundary conditions. In particular, the potential 
V is real. (Absorption may be described by a complex U, but causality then requires dispersion; the necessary 
generalization of (IT) will not be discussed here. 



The complex QNM frequencies are often directly observable: for example the central frequency and width of an 
optical line observed from a laser cavity, or the rates of repetition and decay of a gravitational-wave signal that may 
within the next decade be detected by instruments such as LIGO [0. By way of orientation. Figure la shows a 
square well (solid line, ignore broken line), and Figure lb shows, on the complex a;-plane, the distribution of NM and 
QNM frequencies (crosses and circle, ignore triangle). Several features may be noticed, (a) Bound states or NMs |^ 
are shown in the upper half-plane; in this example there are three, (b) QNMs have Ithuj < 0, and are in the lower 
half-plane. Provided Re a; ^ 0, they occur in pairs: cj_„ ~ ^^ni ^ is readily shown by conjugating the defining 
equation and boundary conditions, (c) There may be QNMs, not paired, which have Rew = 0; in this example there 
is one. These zero modes pT| ] will be of particular importance below. 

Even though the QNM eigenfunctions are not square-integrable and do not form a Hilbert space (at least not in 
the conventional sense), they are nevertheless useful for analyzing outgoing waves. As stated above, QNM frequencies 
are directly observable, and as illustrated by the example in Figure 1, the spectrum is typically much richer than that 
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for NMs. In fact, the spectrum is so rich that, under some broad conditions, namely that V{x) vanishes identically 
outside an interval [—a, a] and has a singularity at a; = ±a, the QNMs together with any possible NMs are complete 
in at least part of spacetime |0 , so that the wave signal observed at a point x can be represented as 

= ^a„0„(x)e-^'^"* , (1.11) 



where the sum is over all QNMs (including zero modes) and NMs |13|. In cases where the QNMs are not complete, it 
may still be possible to characterize the remainder, which could be, for example, a power law in t ||l^,|l^. Moreover, 
when the QNMs are complete, one can set up a formalism that almost completely parallels the case of NMs 
under a bilinear map which is the analog of the usual inner product (but which is linear in both vectors rather than 
linear in one and conjugate linear in the other), the Hamiltonian turns out to symmetric (the analog of self-adjoint). 
This allows much of the well-known mathematical formalism developed for NMs to be transcribed. The mathematical 
structure can also be written elegantly using a biorthogonal formalism jl^] . One can even second-quantize using these 
QNMs as a basis (e.g., to discuss photon creation and annihilation in an optical cavity) jl^. These developments have 
been reviewed in Ref. ||l9[]. An important point of this paper is that much of the mathematical structure is preserved 



under SUSY (Section VIl) 



It is straightforward and indeed natural to include total-transmission modes (TTMs) in the discussion as well; these 
are defined by the incoming-wave condition at one extreme and the outgoing-wave condition at the other. The TTMs 
propagating from the left (TTMl) and those propagating from the right (TTMr) satisfy 

(t){x) (X 6=^'""^ , ^ OO , (1.12) 

with the + (-) sign for a TTMl (TTMr). 

Given the striking parallel between NMs and QNMs, and especially since the QNM spectrum is often richer than 
the NM spectrum, it is logical to ask whether the entire concept of SUSY can be generalized to QNMs: namely, given 
a real potential V, can one find a real partner potential V such that the outgoing waves for the two have the same 



set of complex eigenfrequencies (or the same set apart from one state) and are related by (1.5)? The same question 
can be asked of TTMs. This paper sets out to answer this question, which turns out not to be more difficult, even 
though one has to match a set of complex rather than real frequencies. Indeed, when NMs and QNMs are studied 
together, the relationship between the spectra of H and H becomes clearer. 

Relations between such pairs of open systems have been discussed in terms of the transmission amplitude T{uj) and 
the reflection amplitude R{io) etc., for solutions at definite frequencies. However, an operator SUSY transformation 



such as (1.5) is more powerful: given a (l}{x,t) that solves the dynamics of H, one can immediately generate the 
corresponding solution (f>{x, t) = A(l){x, t) which solves the dynamics of iJ, without having to first project (j) onto 
frequency eigenfunctions. 

SUSY is closely related to the question of inversion: given the complete set of complex QNM frequencies, can V{x) 
be determined uniquely? Obviously, if isospectral SUSY partners exist, then unique inversion would not be possible. 
The converse need not be the case, since there could be isospectral systems not related by a SUSY transformation 



such as (L5). It may be well to recall the classic results for NMs [^: for a finite interval, say [0, a], if the (real) 
spectrum is given for the two dynamical systems defined respectively by (a) (/)(0) — (f>{a) — (the nodal system), and 
(b) 4>'(0) = 4>{a) = (the antinodal system), then the potential can be uniquely determined. (Equivalently, one can 
extend the problem to [—a, a] and seek a symmetric potential V{x) given the spectrum fo r (a) t he even-parity sector 
and (b) the odd-parity sector.) Effective numerical algorithms have recently been given |2l|,^. Since the QNMs of 
open systems are complex and carry twice the amount of information, it is tempting to speculate that the frequencies 
for one set of boundary conditions (e.g., outgoing waves at both ends of a finite interval) may be enough for unique 
inversion. There is some numerical evidence to support this conjecture, at least in some cases p^j2^ . 

The possibility of inversion from one spectrum is intriguing, since physically all the information one has is the 
time-dependent signal and hence in principle the set of complex frequencies under one set of boundary conditions 
psf . For example, without wishing to belittle the technical difficulties involved, one may imagine that the signal 
received by LIGO could yield the complex eigenvalues and from these one could determine V{x), which describes the 
background curvature. This could then, in principle, provide a novel astrophysical tool — using gravitational waves to 
probe the intervening spacetime between the source and the observer. It has been suggested that the QNM spectrum 
may be a definitive signature of black holes [ p6[ . Moreover, a toy model has been given for the corresponding forward 
problem of perturbing a black hole: if the hole is "dirty", i.e., perturbed by surrounding matter, how would this be 
revealed in the shifts of the complex frequencies p7|? The inverse of this perturbative problem would in principle 
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allow the nature of the astrophysical environment around a black hole to be revealed through the shifts in the QNM 
spectrum. 

With these motivations, in this paper we study SUSY for QNMs and TTMs. The formalism is given in Section ||, 
and is in fact almost the same as that for NMs |28]. SUSY transformations can be classified into three discrete types 



plus one continuous type. Examples for square wells and barriers are given in Section III. It is found that the familiar 



transformations not only establish an equivalence between the NM spectra of partner potentials, but preserve the 
QNM spectrum as well. Moreover, by considering NMs and QNMs together, a clearer picture emerges. The reflection 
and transmission amplitudes for SUSY partners are related in a simple way, as is shown in Section IV . Section ^ deals 
with the Poschl-Teller (PT) potential, which has a number of interesting features: (a) for a range of parameters, the 
FT potential has an infinite number of zero modes, and many of these allow corresponding SUSY transformations, 
leading to infinitely many essentially isospectral partners; (b) two of these transformations lead to partner potentials 
also of the PT type, but with different amplitudes; (c) thus the partner can again be transformed to yet other PT 
systems, ad infinitum; and (d) a subset of PT potentials are SUSY-equivalent to the free field, and therefore have total 
transmission at all positive energies. In Section Vl, the Regge- Wheeler (RW) and Zerilli potentials, which respectively 
describe the axial and polar sectors of linearized perturbations of a Schwarzschild spacetime, are shown to be SUSY 
partners. This result is not novel, b ut is here cast into the general framework of SUSY. The orthonormality of QNMs 
under SUSY is presented in Section VII, in which SUSY is generalized to the two-component formalism appropriate 
to open systems. SUSY can also lead to doubled QNMs and TTMs, in the sense that two such modes occur at the 
same frequency. This phenomenon, which has no counterpart in conservative systems, is analysed and illustrated in 
Section VIII. Finally, a discussion is given in Section IX. 



II. FORMALISM 



It is readily shown, without reference to any boundary conditions and hence with equal validity for NMs, QNMs 
and TTM s, th at H and H have the same spectrum (with the possible exception of one state, to be discussed below), 
provided (1.7) holds. It is readily shown that the two potentials can be written as 



v{x) = w{x)^ - w'{x) + 

V{x) = W{xf + W'{x) + 



(2.1) 



with W{x) (called the SUSY potential) as in (1.6) and for some constant fi^. Moreover, the Hamiltonians can be 
represented as 



H = AA^ + 



(2.2) 



It follows that if (j){x) is an eigenfunction of H, then the partner function (j){x) = A(j)[x) (provided it does not vanish) 
is an eigenfunction of H with the same eigenvalue. The two partner systems can be put into one linear space by 
introducing Pauli spinors, with H and H associated with 1 ± CTz and A, associated with (j±. This transcription of 
the usual formalism remains valid even when the eigenvalues are complex. We defer all discussions of normalization 



to Section VII 



Upon reversing the sign of PF, (a) V and V are interchanged (see (2.1)), and (b) A and —A^ are interchanged (see 
(1.6)); thus the inverse mapping from H back to H is (up to a sign) achieved by . 
We may regard ( ^.ij ) as a Riccati equation for W in terms of the given V . As V{x) is assumed to vanish for 



±oo, we have W — —VL as x ±cxd. The first-order Riccati equation can satisfy two boundary conditions only 



at special values of ; this condition becomes familiar if we define the SUSY generator $(a;) by 



W{x) = - 



<^'{x) 
$(.t) 



Then (2.1) implies that $ must be an eigenfunction of H with eigenvalue fi^: 

H^{x) = n^^{x) . 



(2.3) 



(2.4) 



The various SUSY transformations are then related, in a one-to-one manner, to solutions of (2.4). Before proceeding 
further, we should note that even for the QNM case with complex frequencies, since V is required to be real, W and 
hence ■!> and fi^ also have to be real (in the case of $ up to an irrelevant overall phase) . 
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The task is therefore to analyze the boundary conditions. First, suppose fl^ > 0, so that f2 is real. Then outside 
the support of V, <I> is oscillatory: either complex (e.g. e*^^), inadmissible since it leads to a complex W; or real (e.g. 
sin rice), inadmissible since its nodes lead to singularities in W. Thus, fi^ < 0, and we denote K = \ The somewhat 
special marginal case = is dealt with in Appendix |B 6| . 

At each spatial extreme {x — *■ ±oo), the solution can be either (a) purely decreasing (<i> (x e"^'^'), denoted as D; 
(b) purely increasing ($ cx e^'^^l), denoted as I; or (c) a mixture of the two terms, i.e., 

$(x) = ce^l^l +de--^l^l , (2.5) 

with c,d ^ 0, to be denoted as M. Considering both extremes, the possible solutions for $ can be classified into the 
following types, in obvious notation: 

Typel =DD , 
Type 2 -II , 
Type 3 = DI, ID , 

Type 4 = MM, MD, MI, DM, IM . (2.6) 

The first three types can only occur at discrete values of il, and can be characterized as follows. A solution decreasing 
at both ends is an NM (DD = NM); we associate these with waves going as e*^''^', with Q — iK. A solution which 
is purely increasing at both ends is a QNM (II — QNM); we again associate these with waves going as e*^'^', but 
now with n = —iK. A DI solution goes as e^^ at both extremes {K > 0). We can assoiciate this with either (a) the 
time-dependent solution e*^^e~'^* with — —iK, namely a TTMl, or (b) the time-dependent solution e~*^ x^-ifi t 
with fl' = iK, namely a TTMr. Purely as a matter of convention, we shall only discuss states in the lower half-plane 
(in this case the TTMl), always noting that there is a TTM propagating in the opposite direction at the opposite 
value of the frequency which is not shown. Thus, with this convention, DI = TTMl; similarly ID = TTMr. This 
accords with the convention in discussions of black-hole modes. We shall in the rest of this section pay particular 
attention to the disappearance and appearance of the four classes of solutions (NM, QNM, TTMl, TTMr) under 
SUSY. 

In contrast, Type 4 solutions are allowed for continuous il. There is an additional difference. For Types 1 to 3, the 
logarithmic derivative outside / is exactly ±X, so W = 0. Then, since V vanishes outside /, so does V. However, 
for Type 4, the logarithmic derivative is given by (e.g., for x —f +oo) 

W^-K+^e-^''- + --- , (2.7) 
c 

so that V would have an exponential tail, taking us outside the simplified theory presented in this section. 

We should check immediately that outgoing/incoming boundary conditions are preserved under SUSY. To do so, 
let (j) be, for example, an outgoing wave, (j){x) = Ce*"^ as a; ^ oo; then operating with A, we find 

4>ix) = [dx + W{x)] (b{x) 

= {iuj + W+)Ce"^'= , (2.8) 



where we have introduced the constants 



W{x^±oo) = W± , (2.9) 



which, depending on which type is being discussed, would equal ±K. Thus </> and always satisfy the same type of 
boundary conditions. 

Two properties of the generator <& ought to be mentioned. First, it is annihilated by the operator A: A^ = 0; this 



follows trivially from ( |1.6| ) and ( |2.3| ). Second, from (2.1) we see that the partner systems are related by reversing the 
sign of W, so in view of (2.3), the corresponding generator for the inverse transformation from H to H is ^ — <I>~^; 
this function is guaranteed to be an eigenfunction of H also with the eigenvalue fl^. (Despite the notation, $ is not 
the SUSY partner of $: $ ^ A^.) Because of the reciprocal relation, the boundary conditions for $ are obvious: 
under $ i-^ l>, D i-^ I, 1 1-^ D, and M i-^ D. 

Thus, under SUSY, one state disappears and an extra state $ appears, with boundary conditions (e.g., DD, II) 
readily deduced from the above description. In view of the association of these boundary conditions with the four 
types of solutions (e.g., DD = NM, II = QNM), we can summarize the changes in the number of these solutions as 
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in Table 1. In this table, the entries ±1 in the "NM" column refer to transitions 0^1 only, because there can be no 



doubled bound states (cf. at the end of the Introduction) in 1-d systems; see also Appendix B 5 



In particular, if $ is itself an NM or QNM (Type 1 or 2), then in the partner system one state at f2 = ±iK 
disappears and another state at — f2 — ^iK appears instead. In these cases, the systems, differing by one state whose 
frequency is reversed, are said to be essentially isospectral, and the mapping A between them constitutes a good SUSY 
transformation Hi. Where $ is not an NM or QNM of H, and l> is not an NM or QNM of H (Type 3), then the 
two systems are strictly isospectral (for both the real NM and the complex QNM spectra), corresponding to broken 
SUSY |§. 

Some of these properties are of course well known; but the more systematic perspective is possible only if attention 
is paid to the QNMs (and TTMs) as well. These general remarks will next be illustrated with examples. 



III. SQUARE WELLS AND BARRIERS 

In this section we illustrate the different types of SUSY transformations using square wells and barriers for V] this 
has the advantage that all different cases can be solved analytically. All examples are presented graphically with a 
uniform format. The top diagrams show the original potential V (solid line) and the partner potential V (broken 
line). The bottom diagrams show the complex w-plane. The NMs and QNMs common to both systems are shown 
with crosses. The mode corresponding to $ ($) is shown with a circle (triangle); these are present in one system but 
not in the other. In each case, we have verified the lowest few QNMs of H numerically. 



A. Type 1: $ is a normal mode 

The Type 1 solution, denoted as DD, corresponds to $ itself being a bound state or NM. This occurs only if V 
is sufhciently attractive. No matter how many bound states there are, only the ground state is acceptable as the 
generator <&, since the others would have nodes leading to a singular W — exactly as in the familiar discussion of 
SUSY between NM spectra. In this case, the two spectra (both NMs and QNMs) are identical except that: (a) H has 
one fewer NM, since the ground state in H , namely 0o — is annihilated by A\ (b) H has one extra QNM, since the 
corresponding function $ — <^^^ is a solution with eigenvalue Vl^ but which is increasing (i.e., is a purely outgoing 
wave) for x — > ±cxd. The systems are essentially isospectral. 

As an example, take y to be a square well (Figure 1): 

V{x)^Voe[a-\x\) , (3.1) 

where 9 is the unit step function, and with parameters say Vq = —20, a = 1. The even wavefunctions, for example, 
are given by 

(l){x) — Acoahax , (3.2) 

for |a;| < a, where 

- a2 _^ Vb = . (3.3) 

The logarithmic derivative at a: = a must match onto an outgoing wave: 

a tanh aa — iuj , (3.4) 

where the case of NMs is included if we take Re w = 0, Imw > 0. The odd sector can be treated similarly. 

For this choice of parameters, there are 3 bound states, at w = 4.28i, 3.68i, 2.47i. Taking the lowest of these as the 
auxiliary function $ (so that f2 = -\-iK with K = 4.28), we find for |a;| < a 

(^{x) = Acosqx , 
W{x) = gtanqx , 

V{x) = Fo + 2q^ sec^ qx , (3.5) 

where q = \/—Vo — is real; with the present parameters q = 1.28. Figure la shows the two partner potentials and 
Figure lb compares the spectra. The NM in H but not H is shown by the circle a.t ui — iK; the QNM in H but not 
H is shown by the triangle at w = —iK. 
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B. Type 2: $ is a quasinormal mode 



The Type 2 solution, denoted as II, corresponds to $ itself being a QNM. Since fl^ has to be real, the frequency is 
purely imaginary and $ is a zero mode. We note two features. First, not all potentials will have a zero mode. Second, 
in distinct contrast to the NM case, there may be several nodeless QNM zero modes, all of which are candidates for 
the generator 4>. 

To illustrate the first property, consider the example of a square barrier, i.e., (3.1) with Vq > 0. Again we take 
a = 1, and allow Vq to vary. The formulas are similar to (3.1)-( [3^ ) with the sign of Vq reversed, and will not be 
shown. In the even sector there are two zero modes for small values of Vq; as Vq increases, these move closer together 
and merge eX Vq = Vc = 0.291 pl| ; for higher barriers there are no zero modes. There are no zero modes in the 
odd sector, and in any event an odd state is not eligible as a generator (because there is a node). To be definite, 
we take Vq — 0.16, in which case the zero modes are at w = — O.lSli, — 2.500i, with a — 0.242, 2.506 in ( |3.2D . These 
wavef unctions, being cosh ax within [—a, a] and a real exponential for \x\ > a, clearly have no nodes. 

This takes us to the second property alluded to above, namely, that the existence of more than one nodeless 
QNM is by no means exceptional; some general remarks are given in Appendix Multiple nodeless QNMs allow 
different SUSY transformations. Continuing with the example, we show separately the case where the state at 
n = wi = — O.lSli is chosen to be $ (Figure 2) and the case where the state at uj2 = —2.500? is chosen (Figure 3). 
The two partner potentials V are different. In each case, one QNM has disappeared and an extra NM has emerged. 
The systems are again essentially isospectral. 



C. Type 3: $ is a total-transmission mode 



To be definite, we take the Type 3 solution to be DI: $ is decreasing as x — oo and increasing as x — )■ cxd, i.e., a 
TTMl according to the convention adopted. A simple example can be provided by a multi-step square barrier: 

f Vi if \x\ < b 

V{x) ^ I Vo if 5 < |x| < a . (3.6) 
[ if |x| > a 

We take b = 0.1, a = 1.0, Vi = -10.0, Vq = 1.0. A TTMl is found for n = -0.990i, and it can be used to generate a 
strictly isospectral SUSY partner: all the NMs and QNMs are preserved. Figure 4a shows the partner potential, and 
Figure 4b shows the QNM distribution. Figure 5a shows the generator, while Figures 5b and 5c respectively show the 
NM (pi and the QNM (j)2 in the partner system V, with eigenvalues oji = 0.498z and lu2 = — 1.570i. (Because in this 
example V is symmetric, there is a TTMr at the same fi. However, that mode would generate a different V, in fact 
the parity-partner of the one considered here.) 



D. Type 4: Continuous transformations 



As an example, consider the same square b arri er as in Section IIIB , i.e., Vq — 0.16, a = 1.0. There are many 
sub-categories for Type 4 transformations (cf. ( |2.6| )). We choose symmetric boundary conditions, i.e., (2_^) with the 
same values of c and d for both x — > ±oo. Within the range |a;| < a, is again given by (3.5). However, for > a, 
there is an exponential tail: 



V{\x\>a) = 2K'^ 



1 - 



^1 - de 



-K[. 



ce 



K[. 



de-K\^\ 



(3.7) 



Turning to the mode structure, it turns out that $ is "universal" in that it is an NM, a QNM, and a TTMl/r at the 
same time, because outgoing and incoming waves coincide in the partner at w = —iK. (This somewhat surprising 
phenomenon is explained in Appendix ^.) As in other cases, modes at eigenvalues 7^ —K^ are not affected by 
SUSY. Figure 6 shows an example, where we have chosen K = 3.0, implying d/c ~ —0.829. 
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IV. REFLECTION AND TRANSMISSION AMPLITUDES 



The relation between the reflection and transmission amplitudes of SUSY partners is known in general [pf , and has 
also been used in the specific context of gravitational waves propagating on a black-hole background H|T Here we 
record the results for completeness. 

The reflection and transmission amplitudes (for a wave incoming from the left) are deflned by 

for X — > — oo a nd a; — > oo respectively. Now applying A = dx + W and recognizing that W{x) for large x equals finite 
constants (2^), the corresponding wavcfunction in the partner system is 

respectively for x — *■ =poo. Normalizing the incoming wave on the left to unity, we find the corresponding amplitudes 
in the SUSY partner to be 

f,(„) = iiil^r.M . (43) 

The constants W+ and W- are real and always equal in magnitude. For real u; it then follows that |i?L('^)P = 
|i?L(w)p, |Tl(ijj)P = |Tl(cli)P. In the case of Type 3 SUSY transformations, W- — W+ and one has the stronger 
condition Ti^{uj) — Ti^{uj). 

Similar relations can be obtained for the amplitudes and defined for a wave incoming from the right; it 
suffices to change W± — > and uj —u) in the prefactors above. In fact, one has the relations 

Ti,{uj)Ri,{-u}) 

TrM = TlH . (4.4) 

Although the transmission and reflection amplitudes for pairs of related open systems are useful, they relate to 
states of dcflnite frequency, and are less general than the operator SUSY transformation. 



V. POSCHL-TELLER POTENTIALS 



A. General properties 

We next illustrate with two examples related to the PT potential |^ : 

V{x) =Vsech^{x/b) , (5.1) 

which has an exponential tail: V{x) ^ Ve^^'^'/''. This means that the theory of Section || is not strictly sufficient 
(cf. Appendix^), but for simplicity we shall gloss over the complications as much as possible, only pointing them out 
where appropriate. The PT potential admits analytic solutions, and because of its exponential tail is sometimes used 
as a proxy to illustrate some properties of the effective potentials that arise in the study of gravitational waves; the 
latter also have exponential tails for a; — > — oo, where x is the tortoise coordinate and — oo corresponds to the event 
horizon. (See Section ^ for details.) The PT potential is especially interesting from the point of SUSY, in at least 
three ways. 

(a) There is a class of SUSY transformations (Types 1 and 2) which map one PT potential to another of the same 
width {V ^V,b^b). 

(b) There are Type 4 transformations which map the free field (V = 0) to a PT potential with b'^V = -2. 

(c) By again transforming this potential as in (a), one obtains a string of other PT potentials with b^V = —£(£+1), 
with £ being an integer. 
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These results are interesting because (a) implies that PT potentials are shape-invariant under SUSY, while (b) and 
(c) show that certain PT potentials are SUSY-equivalent to a free field, and therefore have total transmission for 
all positive energies. Attractive PT potentials have a finite number of bound states, while repulsive PT potentials 
obviously have none; thus the existing SUSY literature, which focuses on NMs, has not given a systematic account of 
these properties. 

To pave the way for discussing SUSY, we first summarize the known results for the PT potential. For any uj, the 
wavefunction that is outgoing at a: — * —cxd reads 



'X 2Fi{^+q~iujb, ^—q—iujb;l—iujb;£^) , 

where 2F1 is the hypergeometric function, and the variable 1^ is related to x by 

1 



1 



-2x/b 



The parameter 



q=Jl-b^V 



will be especially important below. For generic q, these solutions are outgoing at the following values of u>: 



'-{n+'^±q) 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



for n = 0, 1, 2, . . ., which thus are the (Q)NM frequencies of the PT potential. The states have parity (—1)". However, 
ii q = 5 + £ (^ = 0, 1, 2, . . .), the modes are restricted to iujb — ~£, —£+ 1, . . . ,t. apparently, for half-integer q the two 
strings of modes (5.5) annihilate each other where they coincide. (Moreover, for these q the remaining modes are of 
a special type, coined "universal" in Section HID; see further Appendix ^) 

From (5.4), it is clear that the value b^V is important, and in fact the situation can be classified as follows. 



-K1 



This case will not be discussed 



(a) If J < b V (so that q is imaginary), the QNMs are paired with lo^ 
here, because there are no real SUSY generators. 

(b) If < b^V < i (so that q is real and the potential is repulsive), these QNMs are all zero modes lying on the 
imaginary axis. This case is interesting because many of these QNMs can be used as the auxiliary function $. 

(c) If b^V = 0, one obtains the free field, whose SUSY partners will be discussed in Section VE . 



(d) li b^V < 0, one obtains an attractive PT potential, and it will be seen that these could be SUSY partners of 
repulsive PT potentials. 

(e) If b^V = —£{£ + 1), the attractive PT potential is SUSY-equivalent to a free field. 

These cases will be discussed in detail below. In fact, some of them could be anticipated from the eigenvalue formula 
(pTq), which has the important property: 



for integers n and n iff 



Likewise 



for integers n and h iff 



q = q± {n — n) 

(9) = (9) 
q = —q ± (n — n) 



(5.6) 
(5.7) 
(5.8) 
(5.9) 



Secondly, although (5.4) implicitly defines q >0 (where it is real), negative values can be admittted as well, and 

w+(g)=<(-g) ; (5.10) 



these are just different labelings of the same state. So by allowing negative values of q, ( ^.8D -(5.E) could be subsumed 
under ( |5.6D -(5.7). The equivalence in spectrum implied by (5.6) and (5.8) (which is not strict isospectrality because 
n and h have to be non-negative) already suggests a SUSY relationship, which will be presented in the following. 
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B. Results for truncated potential 



Because of the subtleties when the potential has a tail, we first consider a slightly altered problem — the PT 
potential truncated at x = \a\: 

Vix) ^Vsech^{x/b)e{a~\x\) . (5.11) 

We take V — 3/16, b — 1, and a ~ 2 (solid line in Figure 7a). The QNMs can be found by requiring the logarithmic 
derivative at a: = ±a to be ±ia;; the resulting frequencies are shown in Figure 7b (crosses and circle only). The 
lowest QNMs are approximately the same as in the untruncated case, though the higher QNMs are altered beyond 
recognition [^3| , while many new QNMs appear because of the reflections a.t x = ±a. The important point is that 
there are still two zero modes, in this instance a,t uj — — 0.224i, — l.SOlz. Either of these can be used as $; to be 
definite, we take the first. The SUSY transformation is thus of Type 2. The partner potential is shown in Figure 7a 
(broken line). We note two features: (a) V vanishes outside [—a, a], and (b) V in this case is attractive. The overall 
spectrum is shown in Figure 7b. Notice, as before, the appearance of an extra NM accompanying the disappearance 
of the QNM. 



C. Results for full potential 

We now return to the full PT potential, the eigenvalues for which are given by (^). For < |, all the QNMs 
ui^ are zero modes. For states with odd n, the node at a; = prevents their use as the generator <1>. If V > 0, all 
the even-n eigenf unctions are nodeless (Appendix ^) and therefore eligible as the generator, which we denote as 
This then leads to an infinite number of SUSY potentials W^, and correspondingly an infinite number of partner 
potentials Vj^. Since we start with a QNM, these refer to Type 2 SUSY. For attractive potentials, the construction 
goes through at least for n = 0: <i>Q is the ground state and hence certainly nodeless, and $o is readily checked to 
have the same property. Some arithmetic leads to the following explicit expression: 

yn{x)= [V + yl±-2B±C±G2(x)+2(B±)2Gi(a;)2sech2(a;/6)+4B±Gi(a;)tanh(a;/6)]sech2(x/6) , (5.12) 

where 

Fn.pi^) = 2Fi{^+q-iuj^b+p,^-q-iuj:^b+p;l-iuj^b+p;£,) , 

A±-(-2n-lT2g)fe-^ , 
n{n =F 2q)b~'^ 



-2n+l ± 2q 



. (n-l)(n-lT2g) 
^" ~ 3-2n±2q ' ^^'^^^ 

Figure 8 shows the PT potential V and some of its SUSY partners obtained in this manner. 



D. The case n = 



We omit a detailed discussion of all the SUSY partners, but those generated with the n = states as the generator 
deserve special attention. From ( 5.13| ), we see that — 0, and consequently Vj^(x) are again of the PT form, with 
the same width b but different amplitudes, given by 



b^V fe^V - 1 T 2g = fo^V - 1 T Vl - 'ib^V 
This can be expressed more succinctly in terms of q: 

q^\l±q\ . 



(5.14) 



(5.15) 



The equivalence in spectrum can now be readily verified from (5.7) and (5.9). In fact, we can use the freedom in the 
labeling convention provided by (5.10) to rewrite (5.15) as q — q± 1. 
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Since the partner V{x) is another PT potential, we can apply these transformations again and again. The result is 
a chain: 

• • • —m+q <->•••<-> — 1+q q 

1+q ^ • • • ^ m+q ^ • • • . (5.16) 

This shows that the repeated application of SUSY in this manner forms a group isomorphic to the integers. The PT 
potentials at both extremes of ( ^.16D become increasingly attractive, with more and more NMs. 

The PT potential is self-replicating in the sense that (some of) its SUSY partners are the same potential with a 
different amplitude. In Appendix O we show that conversely this condition leads uniquely to the PT potential. 



E. SUSY partners of the free field 



The free field V{x) = is a special case of the PT potential, corresponding to q = i. Apply a Type 4 transformation 
by using the generator 



A simple calculation then gives 



$(a;) = coshi^x 



V{x) = -2K'^sech^ Kx 



which is of the form (5.1). We see that an attractive PT potential is a SUSY partner of the free field if 



(5.17) 



(5.18) 



(5.19) 



or equivalently q — ^- This potential has one NM, which is just $ = 

For half-integer q the QNM oJq is missing, see below (5.5). However, its limit as g — > ^ + i, being the even solution 
(cf. Section HID) at eigenvalue b^il^ = —{£ + 1)^, still generates a SUSY transform which can be verified to increase 
q by one. Successively applying these (Type 4) transformations one obtains increasingly attractive PT potentials, 
described by b^V — —2, —6, — 12, . . . , —£{£+!), g=|,|,|,..., £+^, ■ ■ ., with the number of NMs being respectively 
1,2,3,...,-^,.... 

These PT potentials are interesting because, being SUSY-equivalent to the free field, they obviously have total trans- 
mission for all complex frequencies. This property has recently been discussed in the context of electron transmission 
through semiconductor devices pil, but here it falls naturally into the SUSY framework. 



F. Coalescence of QNMs 

The eigenvalue formula (|]J) implies that QNMs coalesce in two ways as the system parameters are tuned; these 
are quite distinct scenarios, and it is well to emphasize their difference. 

First, as g ^ 0, the QNMs ui^il) ^^"^ ^nil) nierge. At the point of coalescence, they become a double pole; this 
is related to an interesting Jordan-block structure in the evolution operator Beyond this point, the two QNMs 
separate in the real direction (i.e., q becomes imaginary). The transition is exactly the same as a harmonic oscillator 



going through critical damping. This case will be further discussed in Section VIII, The same occurs for q ^ n. 

In contrast, as g — > i, the QNMs ^+((7) and ^,7+1(9) rnerge. However, when they coalesce, the two poles disappear. 
With all poles eliminated, the dynamics become essentially trivial, with total transmission. The same occurs for 



VI. REGGE-WHEELER AND ZERILLI POTENTIALS 



An application of particular interest concerns linearized gravitational waves propagating on a Schwarzschild back- 
ground. In each angular- momentum sector I > 2 and for each spin s (for gravitational waves, s — 2), these waves can 
be described by a scalar field cj) of the radial variable, satisfying a KGE with the potential V describing the nontrivial 
background [Q. However, the axial sector is described by the RW equation, which is (1.1) with the potential 
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V{x) = 21 



2m\ x,(n+l) 



3m 



while the polar sector is described by the Zerilli equation, which is ( |l . 1| ) with the potential 



, , , 2m 
Vix) = 1 - — 



2n^(n+l)a;^ + Qn^mxl + ISnrn^x, + 18rn^ 



x^XnXi, + 3m)2 

Here m is the black-hole mass, and x is the tortoise coordinate, related to the circumferential radius by 



1 



Note that a;* e (2m, +oo) maps to a; G (— oo, +oo). The constant n is related to the angular momentum / by 

n^-{l-l){l+2) . 



(6.1) 



(6.2) 



(6.3) 



(6.4) 



It has been noticed that the two equations have the same QNM spectrum. For example, Figure 9 shows the 
distribution of QNM frequencies for Z = s = 2 for both equations. This can be understood because the QNMs of the 
two systems are related by 



n(n+l) 3m(a;* — 2m) 



— ( 

dx 



3m x1{nx<, + 3m) 

which has been referred to as "intertwining" |Q . Eq. ( |6.5[ ) is exactly a SUSY transformation, if we identify 

2m) 



W{x) 



n(n+l) ^ 3m(x* 



3m xj(nx* + 3m) 

The corresponding generator $ must be an eigenfunction of the RW equation with eigenvalue , where 

.n{n+\) 



3m 



(6.5) 



(6.6) 



(6.7) 



This is exactly the so-called algebraically special frequency |37|-|39t. With the same eigenvalue, $ — <I>~^ should be 
an eigenfunction of the Zerilli equation. Eq. (^^) shows that W+ = ifi; hence, $ is incoming from the right. Since 

H; this can easily be verified using the Leaver series solution |40|), this 
as has sometimes been suggested pi]. However, W- = ifl does not prove 



is not special for a; — > oo (cf. Appendix 
already establishes that $ is not a QNM — 



conclusively that it is outgoing to the left (cf . Appendix B 1 ) , as implied by a conflicting suggestion that <I> may be a 
TTMr In fact it can be shown that $(x^— oo) is not outgoing, so that $ is not a TTMr and (3.5;) is a Type 4 
transform, being of category (b2) on the left in the terminology of Appendix B 3. However, these and other aspects 
will be reported separately. 



VII. ORTHONORMALITY 

A. Orthonormality for NMs 

In the familiar discussion of SUSY for NMs, orthonormality can be preserved. There are two issues: (a) orthogonality 
is preserved because the transformed NMs are eigenvectors of the self-adjoint operator H; and (b) normalization is 
preserved if the transformation is changed to 

(/>ri 1-^ = N„4>n = NnA(j)n , (7.1) 

for each eigenstate n other than the ground state <&, with 

^ {(t>n\(i'n) 2) 
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It is readily shown that N^^ = — r2 , a result that can also be read off as a special case of the derivation below for 



QNMs. Since — is the eigenvalue of the operator A'' A (see (2.2)), the normalized SUSY transformation (7.1 
can also be written in operator form 

^^l = A{A^Ay'^^^ , (7.3) 

valid for any state (j), not just frequency eigenfunctions. In this and similar formulas below, it will be understood that 
the domain is restricted to the subspace orthogonal to on which therefore A'^ A is non-zero. This makes it formally 
easy to verify the preservation of inner products in the mapping from this subspace: 

= {m ■ (7.4) 

However, when operating on a general wavefunction 0, the factor (A^A)^^/^ can only be evaluated by projecting (f> 
onto the eigenfunctions, and scaling each component by Nn. Thus, in practice, the significant result is the evaluation 
of the factor Nn- We now generalize these concepts to QNMs. 



B. Normalization and inner product for QNMs 

It is necessary to first digress and briefly review the concepts of orthogonality and normalization for QNMs. The 
central issue is that with the outgoing-wave boundary condition, H is not self-adjoint in the usual sense, and different 
QNMs are not orthogonal under the usual inner product. Likewise, the norm J \(j)n\'^dx is divergent for a QNM, since 
the wavefunction grows exponentially at spatial inflnity. 

An appropriate normalizing factor for QNMs was first introduced by Zeldovich | |4l[ | , and later generalized and applied 
to other situations ]42| , including models of linearized waves propagating on a background that is a Schwarzschild 
black hole plus a perturbation : 

((/>n,^n) = 2tJ„ / (f>„{x)'^dx 
J —a 

+ i[cj,n{-a)^ +Maf] ■ (7.5) 
This expression goes as (j)^ rather than |(/)„P, and is in general not real. The upper limit of the integral and the first 



surface term can be shifted to any bi > a without affecting the value of (7^); likewise, the lower limit and the second 
surface term can be shifted to any 62 < —a. This expression is the correct normalizing factor in the sense that, e.g., 
under a perturbation V 1-^ V + AV, the complex eigenvalues of the QNMs change by 

((Pn,<Pn) 

Since one no longer has positivity, there is the possibility that (0„, (/)„) = 0. This exceptional case can be separately 



taken care of |3l|] , and some interesting aspects are dealt with in Section VIII , 

To go beyond the normalizing factor and discuss the analog of orthogonality, one has to first regard each state as 
a two-component vector ip = {ip^,tp'^)'^ = {i/j , dtij')'^ , which is most easily motivated by noticing that the dynamics 
requires two sets of initial data. In terms of the two-component vector, one can define a bilinear map |]^, |l6| , |l9| 

(1^,0) = J jy [tp\x)(P^{x)+ij\x)(t)^{x)] dx 

+ [V;i(-a)</.i(-a) + V^(a)0i(a)]| . (7.7) 



This will be seen to take the place of the usual inner product. Note that for an eigenfunction, (</)„, (/)„) agrees with 
( [7. 51 ). The dynamics can be written in the first-order form idt4> = T^4>^ with the 2x2 evolution operator 

W = *f.02\. ■ (7-8) 



The important property is that Ti is symmetric under the bilinear map: 
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(7.9) 



in the proof of which the surface terms generated in the integration by parts exactly cancel against those in ( [7.7| ) 
The relation (7.9) is the analog of self-adjointness, and leads to the usual proof that for two eigenvectors, 



(0m, 0n) = 



(7.10) 



whenever ujm ^ Un. Provided that (0„, 0„) 7^ [|l|, one can normalize these eigenfunctions in the usu al wa y, i.e., by 
requiring ( 7.10| ) to be 2uim5mn in general. We shall henceforth call this property orthonormality (and (7.1C) alone as 
orthogonality), with the understanding that reference is made to the bilinear map ( |7.7D rather than the usual inner 
product. It also follows trivially that, provided this orthonormal system is complete (which is the case under fairly 
broad assumptions, see Section |), then the time evolution of a wavefunction tjj^t) is given by 



■4>{t) = a„0„ I 



(7.11) 



where the a„ are obtained by projecting the initial data: 

_ (0„,^(^=O)) 
((/)„, 0„) 



(7.12) 



The preservation of orthonormality under SUSY should therefore be sought in terms of the bilinear map (7.7). 



C. Normalization factor for QNMs under SUSY transformation 



We first present a derivation of orthonormality that does not explicitly require the two-component formalism. With 
orthogonality in the sense of the bilinear map already guaranteed by (7.10), it remains only to compute the normalizing 
factor 



[ipnAn) = 2w„ / [{dx + W)(l)n] dx 
J —a 

+ i [0„(-a)2 + <^„(a)2] . (7.13) 

Integrate by parts to convert {dx4>)'^ to —(9^0)0 plus a surface term, express the second derivative in terms of V—uj"^ by 
means of the eigenvalue equation, and write the potential as = — W + 0?. Then, apart from a term proportional 
to oj^ — il^, the integrand becomes a total derivative dxiWcf)^). Using IU(±a)^ = —51^ and dx(f>n{^ci) = ±ici;„(/)„(±a) 
then leads to 



{4>n, 4>n) 



(7.14) 



Incidentally, the case of NMs on [—a, a] (which corresponds to nodal conditions at the ends of the interval) can be 
recovered by simply dropping all surface terms. 

Since the ratio (7.14) takes the same form as for NMs, and is the eigenvalue of A'' A, we can again write the 
normalized SUSY transformation for each eigenfunction as (7.3). 



D. SUSY transformation for two-component form 

We have so far deliberately avoided the explicit two-component form of the SUSY transformation, since, for an 
eigenfunction one can always express dt(l)n in terms of </>„. However, in order to perform SUSY transformations 
on general outgoing wavefunctions (e.g., when given a time-dependent state, to find its SUSY partner at all times), 
the second component must be considered. 

Since the SUSY transformation must commute with time-evolution and 0^ = dt(t>^, it is clear that both components 
must be transformed in the same way. Thus, the (unnormalized) SUSY transformation on the two-component vector 
is ^ = diag {A, A) , which satisfies 
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(^tA,0) = (V',-4t0) , (7.15) 
where = diag (^'f, ^'1'). In deriving ( 7.15| ), one has to integrate by parts. In the integrand, A turns into because 



A contains a single derivative dx, which reverses sign upon integration by parts. The surface terms are seen to work 
out by using, e.g., (/)^(a) = dt4>^{a) = —dx4>^{a), and the known values of W±. 

Note that A'' A = (H — ^1^)1, .4.4^ = {H — 51^)1, i.e., these products do not relate to the two-component H. 

With ( 7.15| ), it is straightforward to show that the normalized SUSY transformation 



4>^ ^ = A{A^A)-'^/^(ly , (7.16) 
defined on the subspace orthogonal to preserves the bilinear map, in a manner that exactly parallels ([7.4[). In 



the exceptional case of SUSYs that generate a doubled mode (see Section VIII), the subspace has to the exclude two 
states on which H — ft^ vanishes. 

These properties show that two apparently unrelated concepts, namely SUSY and the linear-space structure for 
QNMs in open systems (e.g., the replacement of inner products by the bilinear map), turn out to be consistent. This 



should not be surprising: the bilinear map (7.7) has an intrinsic geometric meaning for all outgoing states, not just 
QNMs H]. 

VIII. DOUBLED MODES 
A. Examples 

In this section, we deal with the possibility that doubled modes may be produced by a SUSY transformation. 
Consider a system H with an NM ^{x) at = iK and a QNM '^(x) at -fl = -iK. Then both ^{x) = $(2;)"^ and 
^{x) = A*(a;) are QNMs of H at -fl. Since there can only be one QNM at any frequency *(x) cx ^{x). The 
proportionality constant can be evaluated by 

$ $(-a) 
*(«) 



$(a) 



= -2irj*(a)$(a) . (8.1) 



Incidentally, the agreement of these two expressions can also be seen without invoking SUSY. First one notes that 
$ and 5", being eigenfunctions of Ti. with distinct eigenvalues, are orthogonal. In the expression for the bilinear 
map, the integral term vanishes because the frequencies are opposite, leaving only the surface terms. Thus one finds 
= $) = i[«'(-a)$(-a) -I- *(a)$(a)]. 

Secondly, consider a perturbation of H which splits the eigenvalues (originally both 51^) of $ and in which case 
their counterparts ^{x) = $(a:)~^ and A'i/{x) must also split into two states. Thus the one state represented by 
(X must be the limit of two states that coalesce when the splitting is switched off, and in some sense to be 
made precise below must count as a doubled state. 

This situation has already been encountered. Consider a PT potential with q = 1. Its modes include (see (^.5[)) the 
NM $ at Wfj" = i/2 and the QNM * at = -i/2, both with eigenvalue = -1/4. Use $ to generate the SUSY 
partner with q ~ 0. Its QNMs $ at lJq and ^ = at uJq both occur at the same eigenvalue —1/4. (Other pairs 
also coalesce, but this particular pair is special in that one of the two states is the reciprocal of the generator.) This 
pair of states form a doubled mode, and we have verified that >P c>c $ (cf. (^^)) indeed holds in this example. 

As another example, consider a system H with a TTMl $(a;) at an eigenvalue fl^. The partner system H must 
have a TTMr ^{x) = ^{x)~^ at the same eigenvalue. But if H should have a TTMp> ^'(a;) also with the eigenvalue 
ri^, then ^{x) = A'i/{x) is a TTMr of H with the eigenvalue as well. As before, we expect these states to be 



proportional, and a calculation analogous to (8.1) shows that in this case 



Z = -2ir2*(-a)$(-a) 

= -2irj*(a)$(a) . (8.2) 

Again, this situation has been encountered. Consider the example in Figures 4 and 5. Since V is symmetric, it is 
guaranteed that for every TTMl $(x) there will also be a TTMr 4'(x) = $(—2;) at the same eigenvalue. With this 



substitution for the condition (8.2) is readily verified for any symmetric potential. 
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B. Double zeros of Wronskians 



In this subsection we make precise what is meant by a doubled mode and in the next subsection we describe the 
SUSY transformation of such doubled modes. We shall do so for QNMs, the case of TTMs differing only by some 
signs. In the original system define outgoing- wave solutions of the wave equation f{uj,x) and (7(0;, x), satisfying 
the boundary conditions 

a;>a) = e*"^ , (8.3) 

where ' ~ dx- (For TTMs, e*"^ on the r.h.s. also for /.) A QNM is a solution that satisfies both the left boundary 
condition (as on /) and the right boundary condition (as on g); in other words, at a QNM frequency, / and g are 
linearly dependent. Thus, the zeros of the Wronskian 

Jq(w) = J(/,5;w) 

= /'(^, x)g{uj, x) - f{uj, x)g'{uj, x) (8.4) 

identify the QNMs. In fact, it is easy to show that the Wronskian is related to the transmission amplitude introduced 
earlier by: 

Now consider the analogous construction in the partner system H, obtained by using an NM $ of iJ as the 
generator. By our convention, $ is associated with a frequency on the positive imaginary axis, and W{±a) = 
-<i>'(±a)/$(±a) = =Firi. The SUSY transformation gives 

f{uj,x) = {dx + W)f{LO,x) 

~g{uj,x)^{dx + W)g{io,x) , (8.6) 

leading to the Wronskian 

Jn{.uj) = {[{dx + W)f]'[{d^ + W)g] 

~ m + W)f][{dx + W)g]'} . (8.7) 

This Wronskian is however unnormalized (as indicated by the subscript), since only Cf and Dg satisfy the normal- 
ization conventions at —a and +a respectively, where 

C = -D=^ . (8.8) 

LO — il 

These normalizing factors have no zeros or poles in the lower half-plane, which is the only region of interest for QNMs. 
Thus the normalized Wronskian in the partner system is 

Jq{io)^CDJ^iLo) 

=^ {lO -n)-^Jn{Lu) . (8.9) 

When (^1^) is written out, some terms cancel by using = 0, the second derivatives can be eliminated by the 
defining equation, and V is expressed in terms of W and fl^. Some arithmetic then leads to 

J^{u) ^ {lu^ - n^) J^{lu) , (8.10) 

Ui^) = ^^^Uio) , (8.11) 



which is also readily derived from (4.3), since (3.5) shows that Jq transforms as under SUSY 



Eq. (8.11) neatly summarizes the properties of Type 1 SUSY transformations (cf. Table 1): the QNMs of H and 



H are the same except that the former has an extra QNM at w = —il. In fact, this formula can be used in the 
upper half-plane as well, where the denominator in (B.ll) indicates that the NM at w = is absent in H. Moreover, 



16 



( [7.14 ) on the change in normalization under SUSY emerges as a simple consequence ol 
is related to the Wronskian by |l9[ 



IC), since the bilinear map 



(/(w„),5(w„)) 



duo 



.12) 



Now specialize to the case where H happens to have a QNM with eigenvalue this means Jq(i^) has a zero at 



51. By (B.ll), Jq has a double zero at this point. This then gives a precise definition of a doubled mode. It 



is clear that upon perturbation, such a double zero would generically split up into two simple zeros, i.e., two QNMs. 
The existence and properties of double (and higher) zeros in the Wronskian for QNMs relate to the iss ue of J ordan 
blocks (JBs), which has been extensively discussed recently [Q, and to which we will return in Section VlllC. 



With TTMs, the situation is entirely analogous, and the double zero refers to a Wronskian defined with boundary 
conditions appropriate to the type of TTMs under discussion. To be specific, consider a SUSY transformation of Type 
3a, which deletes a TTMl and produces a TTMr (by convention we only consider TTMs in the lower half-plane). 
We define Wronskians using the boundary conditions appropriate to the latter, and in exactly the same way arrive at 
( ^.11 ). To interpret the prefactor, we note that a TTMr at ±57 is the same as a TTMl at =fO; thus ( ^.11 ) in this 
case simply states that a TTMl in the lower (i.e., a TTMr in the upper) half-plane is deleted while a TTMr in the 
lower (i.e. a TTMl in the upper) half-plane is produced. 



SUSY transformations of Wronskians are discussed in more detail and generality in Appendix B 5 , in the context 
of potentials with tails. 

The coalescence of two modes is similar to degeneracy (though one of the degrees of freedom is not an eigenstate). 
Usually, a degeneracy is associated with some symmetry. It is therefore interesting to note that in our example of 
doubled TTMs, the TTMl and TTMr are at the same frequency because of parity; under SUSY, they then turn into 
a doubled TTMr in the partner system (even though the latter has no parity invariance). 



C. SUSY partner of a Jordan block 



Returning to the situation described at the beginning of Section [VIII A| , let us examine in some detail the relationship 
between the doubled QNM at —Q, of H and its SUSY pre-image. 

To conform with the notation of Ref. ^l[, denote the state ^ as '^j, where j is a QNM index. Now, using 



the normalization convention (3^), one has j{x) — ^'j(— a)/(— 57, x)e*^' 
(Analogous formulas involving $ and $ follow by the proportionality (S 
a double zero at bj = —57, f(uj,x) satisfies the outgoing condition at x 



implying *j(.t) = «'j(-a)/(-57, 



1) and will not be given.) Since Ju(w) has 
= a not only at uj 



-57, but also to first 

order away from the pole. This makes it plausible that, in the QNM expansion, 9t^/(w, a;)|_o takes the place of the 
"missing" eigenfunction when l> and coincide, which has been confirmed in detail |3|]. One thus defines, for an 
arbitrary a, a pair of functions ^'j,n, where the QNM index j labels the JB and n = 0, 1 is an intra-block index: 



%i(x) = *(-a)a./(a;,x)Lj,e'^'^ + a^.ix) 



(8.13) 



where the second function satisfies 



(J? - 572)%i = -257*j-| 



(8.14) 

Using this, one verifies that a time-dependent outgoing solution is ^j^i{t) = — it^j.o)e*^*, showing that the 

associated momentum reads ^ = — «(— 57\l/j^i -I- ^j,o)- One of the issues to be investigated is the behaviour of the 
term te*^* in i under (inverse) SUSY, since such a time dependence cannot occur for H. 

Next consider the bilinear map for these functions; because the second components are nontrivial, we explicitly 
use the two-component notation. For a double pole one always has (^j,o,*j,o) = 0, cf. (S.12). The undetermined 



constant a in (B.12) can be used to also achieve (\I'j_i, \I'j_i) = 0; this is useful in wavcfunction expansions, but for 
our purposes a can be left arbitrary. 

The above merely applies standard JB theory to the double pole in the _ff-system. Turning now to SUSY transfor- 
mations, one can trivially write 



%i(x) =«'(-a) X d^iid^ + W)f{LU,x)]^ 
= {d, + W)^i-a)d^f{u;,x)\^^e 



.15) 
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Thus, 'i'{—a)du)f{u},x)\_^e^^'^ is the SUSY pre-image of ^^.1(2;). However, the former is not outgoing, since in the 
iJ-system the Wronskian only has a first-order zero at —fl, and consequently f(uj,x) satisfies the outgoing- wave 
condition at x = +a only at — J7, and not to first order away from the eigenvalue. 

The normalization of a JB is determined by one overall factor, which equals the bilinear map between the two basis 
states pi|. In the present case it is evaluated as 



-a 



= -2n , (8.16) 



where in the numerator of the first line we have used the result analogous to ( 8.12| ) for a double pole |^ 



Let us finally consider the reverse transform generated by —A^ [Q. It is easy to establish the following three 
properties, (a) A^'^jfi cx yl^$ — 0. (b) Hence in A'^^j,i{t), the term cx te*^* is annihilated, (c) The remaining term 
in A''^j^i is c^j. The last property is readily seen by observing that 

{H - n^){A'''^j^i) ^ A\H - 1 

= A\-2mjfl) 

= , (8.17) 

so that A'^^j i is an eigenfunction of H with eigenvalue VP' . We have also noted that the SUSY pre-image of ^ j^i 
(i.e., i9cj/|_si) is not an outgoing solution at all. These remarks completely resolve the puzzle related to the prefactor 
t in the time evolution in the iJ-system. 

It is instructive to compare various calculations of the constant of proportionality c above. Defining — A^'^ j^i 
(to avoid double stacked tildes), in the two-component formalism one can calculate 

= (*,,!, [if 

= -2^l(*,-i,#,-o) 

= 4172(*j-,*j) . (8.18) 

One can get the same result from (|8.10| ), for which one has to realize that in Ju(ij-') = (j-o^ — (— il)^) Jq(a^), one overall 
factor (w -(- ri)^ is associated with A'^ J{—VL,x) = and should be removed. Thus, one should evaluate 

(*„ *,) _ / -d^[[u + n)-2 J,(w)] 1 



= ~2VL . (8.19) 

For simple poles, one also knows the sign left undetermined by the bilinear map: 0„ = A'l'(^„ = A^A0„ = (cj^ — ri^)(/)„. 
For JBs, the corresponding calculation is 



n 



= *(-a)9^[AtA/]_^e^"" 

= *(-a)9^[(i?~17^)/]_^e*"'^ 

= -2n^j . (8.20) 

Incidentally, if the eigenvalues uij of the JB and fl^ of the SUSY auxiliary function are not equal, the transformation 
leaves the JB structure unmodified. An example of this situation has already been encountered in the PT potential 
with q = 1, for the QNMs with usj 7^ — i. We omit the details, which are straightforward, as is the generalization to 
higher-order multiple poles. 
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IX. DISCUSSION 



In this paper, we have extended the usual discussion of SUSY as a relation between NMs of partner systems to 
include the QNMs as well; TTMs also come in naturally. By viewing all these together, a more complete picture 



emerges, most conveniently summarized by Table 1 or say (8.11). For example, in the usual discussion for NMs only, 



essentially isospectral systems are described as being identical except for one NM present in one system but absent in 
the other; now we see that (under Types 1 and 2 transformations) when an NM appears (disappears), a corresponding 
QNM disappears (appears). 

QNMs differ from NMs in two regards. First, they have complex frequencies; nevertheless, even with twice as 
many constraints, matching the spectra is not any more difhcult. Second, QNMs need not have an increasing number 
of nodes, and it is often possible to find several nodeless QNMs which generate distinct SUSY transformations — 
whereas the analogous operation for NMs would be restricted to the nodeless ground state. 

The general formalism also provides a framework in which a number of well-known results find a convenient and 
unified expression. These include (Types 1 and 2) SUSY transformations relating different PT potentials of the same 
width, as well as (Type 4) SUSY transformations that take the free field into a nontrivial PT potential. The latter 
property provides one way of understanding the total-transmission property (at all positive energies) of a class of PT 
potentials. The transformation between the RW and Zerilli potentials can also be placed into this framework. 

This wider perspective is gained only because attention is paid to the Klein-Gordon rather than the Schrodinger 
equation, since the concept of outgoing waves has no meaning in an equation that is first-order in time. 

SUSY also preserves the orthonormality of QNMs in the sense of the bilinear map — supporting the view that both 
the extension of SUSY and the generalization of the inner product are sensible and useful. 

Two further important properties are also preserved, (a) If V has a singularity say at x = ±a (e.g., a step), then 
V will have the same type of singularity (but with the opposite sign). This is illustrated by many of the figures, and 
can be seen from ( ^.l[ ) by noticing that the most singular part is W . (b) If V has no tail, then for transformations 
of Types 1-3, V likewise has no tail. These two properties are precisely the necessary conditions for the QNMs (plus 
any possible NMs) to be complete [|l^. Thus, SUSY maps a complete basis to a complete basis, if for Types 1 and 2 
$ i-^ $ is included as well. 

Finally, this work partially answers the question of QNM inversion: does one set of QNM (again plus NM) fre- 
quencies (roughly speaking carrying twice the information of one set of NM frequencies in closed systems) uniquely 
determine V7 In general the answer is negative, for there can be strictly isospectral potentials: see Types 3a and 3b 
in Table 1, and also Figure 5. However, if we consider a half-line problem a; > (say the radial variable in a 3-d 
system), imposing a nodal condition at x = and the outgoing- wave condition at a; > a, can one set of QNMs 



uniquely determine V7 SUSY transformations (1.5) do not directly rule out this scenario — for which there is some 
numerical evidence [|2^ — since these one-sided systems do not feature TTMs generating strictly isopectral partners. 
Moreover, by (pTl|) and (2.3) the nodal condition maps a regular V to V ^ for a; ^ 0+ (generalizing the result 



that SUSY increases the angular momentum by one unit in the hydrogen atom). It would therefore be interesting to 
see if an enlarged class of transformations can address this question. 
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APPENDIX A: TRANSFORMATION BETWEEN THE WAVE AND KLEIN-GORDON EQUATIONS 

In several previous papers |l2| , we have already commented on the relation between the KGE studied in the main 



text and the wave equation (WE) (Al) below. However, this relation exhibits some subtleties in the presence of bound 
states, or when the domain is the full line. Since the main text has shown both of these possibilities to be relevant to 
SUSY, in this appendix the transformation between the two equations will be investigated in detail. 
Consider the WE 

[piz)d^~dMz,t) = . (Al) 
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In mechanics and optics, p is the mass density and the refractive index squared respectively, and thus is positive. 
Violation of p > would also make the time evolution singular and would take one outside the class of hyperbolic 
equations; we the refore exclude this. Then, one can define x by dx/dz = n{z) = \J p{z), 4> = \/nip, and 



upon which <f>[x,t) is seen to satisfy the KGE (1.1). The transformation determines x up to an immaterial additive 
constant, which can for instance be fixed by letting z — Q map to a; = 0. Clearly, if p{z) 1 for \z\ oo then 
V{x) — > for |a;| — > oo (excluding pathologies). 

The WE I— > KGE transform thus is well-behaved at infinity, but singularities in p require caution. While a kink in 
p leads to an admissible S-peak in V, already a step in p is too bad: J V dx = J Vndz = n'/2n^ + J {n'Y /An? dz, with 
the second term diverging if n tends to a step through a sequence of smooth functions. The difficulty stems from the 
fact that both the WE and the KGE become distributionally undefined if p (or V) is more singular than a (5-function. 
The latter is the marginal case, in which the (5-function multiplies a wave function which is not differentiable but still 
continuous. 

Consider the two-component formalism. Since the WE i-^ KGE transform leaves dt invariant, one has cfp' ~ cf?- = 
^/n^^ = n~^/'^ilP'. For the bilinear map this leads to 

(0i>2)kge 



.J — a 



= (^l,^2)wE > (A3) 

where a and (3 are the images of —a and a under the z ^ x map respectively, and where in the surface terms we 
have used the fact that n must be continuous for the transformation to exist. Thus, the WE i— > KGE map leaves the 
linear structure for open systems invariant. 



When investigating the inverse transformation, the form ( A2 ) is not immediately useful since there is no variable z 
to start with. Simple rewriting yields Hq = 0, dz/dx = q{x)^, (j) = qip with p = n? = q^. The boundary condition 
to be used is lim^-^oo 9 = 1- As long as q > 0, one sees that n > and hence p > as well. However, at a zero 
of g, dz/dx suffers a non-integrable singularity and the transformation breaks down. For V > Q (in field theory, in 
the absence of broken symmetry V = m? > 0) this never happens, while if V{x) < for some x, the transformation 
always exists locally but not necessarily globally. 

More precisely, the transformation exists iff H has no bound state. (For the KGE, the main text shows that these 
lead to outgoing solutions which grow exponentially in time, evidently having no WE counterparts.) Namely, the 
breakdown of the transformation means that the bounded lu — wave function has a node, implying that there is an 
integrable eigenstate with lower w^. 

For a KGE defined on the positive half-line (with a nodal condition on the fields at the origin), naturally the case 
in which q{x) has a single node at x = is marginal, and < x < oo can just be mapped onto the entire z-axis. For 
problems defined on the full line, another subtlety shows up for x — s- — oo, where even in the absence of nodes in q the 
transformation is regular only if H has a mode at w = (cf. Appendix [B 6| ), i.e., if q tends to a constant for x —f — oo 
as well. For all other V (which includes all nontrivial nonnegative potentials), q grows linearly for x — *■ — oo, and the 
X-axis maps to a semi-infinite z-interval with a singularity in p at the lower boundary. Of course, by supplementing 
the differential equation for q with different boundary or initial conditions, in this latter case one can also map the 
X-axis onto the negative-z semi-axis, or onto a bounded z-interval. 

APPENDIX B: POTENTIALS WITH TAILS 
1. Preliminaries 

Determining outgoing waves when the potential has a tail is delicate. Namely, whereas the boundary condition 



(1.10) is readily implemented at finite x = ±a whenever these enclose the support of V, for potentials with tails the 
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condition can only be imposed asymptotically. However, if Imw < then lima;^oo 4''{^)/4'{^) = for o.ny solution 
at eigenvalue except the (unique up to a prefactor) small one, so that this criterion does not identify the outgoing 
wave. 

If V vanishes at infinity faster than any exponential, one can impose the boundary condition at some a; = L to 
find an approximate outgoing wave, which will converge to the exact one as L — *■ oo. In practice, for any required 
accuracy one thu s can simply truncate V , so that this case is still essentially the same as the one with no tail. 

If however V{x) ^ e~^^, the above will only work for Imu; > —A/2, thus identifying only QNMs within a band 
below the real axis. By using the mth-order Born approximation for the asymptotic region, one can expand this band 
to imw > -(to+1)A/2 Q. 

For more heavily damped QNMs and for potentials which vanish sub-exponentially (i.e., A = in the above), 
sometimes (semi-)analytic solutions enable one to continue the outgoing wave from the upper into the lower-half 
w-plane, cf. Section refonewave. Fortunately, this is the case for the PT, RW, and Zerilli potentials. In general, 
however, one needs sophisticated numerical techniques such as phase integrals, which involve integrating the KGE 
along carefully chosen paths in the cc-plane citeand. 



2. One outgoing wave 



For general potentials the calculation of outgoing waves thus can present difficulties; in fact already their definition 
needs some care. Therefore define solutions of the KGE f{uj, x) and g{uj, x) such that for Imw > 0, /(w, x) ~ 1 • e~'"^ 
for X — > —oo and g{uj,x) ^ 1 • e*"^ for x — s- oo. This makes / and g analytic in this half-plane, and they can be 
analytically continued to the whole w-plane, with a cut on the negative imaginary axis. Functions oc / and oc g are 
said to be outgoing to the left and right, respectively, at frequency lo; a function is incoming at lo iff it is outgoing at 

— LO. 

As a consequence, one will have g{—uj*,x) = g{uj,x)* and g{u}, x—^oo) ~ 1 • e*'^^ for generic points in the whole 
w-plane. It can happen that the latter relation does not hold at isolated points lo' w here g diverges; we will call 
such frequencies special points and write locally giu) = (lu — cli')^"x('-^)j with x{uj) finite near lu — to' and x(^') 7^ 0: 
and with an integer n > 1, to be called the index of the special point. Near the (possible) branch point lo = the 
behaviour of g can be more complex, and it also can happen that limt^^o di^) = 0. 

From the Born series for g one sees that, if F(x— >oo) ~ e^^^, frequencies lo' = —ikX/2 {k = 1,2,...) can be 
special points: depending on higher-order contributions to V, the divergence in g could cancel. Examples of such 
"miraculous" cancellations include the PT potential, where the special points generically present in /(w,x) given by 
(5.2) a.t uj — ~in/b [n ~ 1,2,.. .) are absent \i q — ^ + I {(. ~ . . .) and n > ^ -I- 1, as follows from standard 
properties of the hypergeometric function p5[ | ; for (, — Q this becomes the free field, which obviously has no special 
poin ts at all. Without proof we mention that such a "miracle" also occurs at the algebraically special frequency 
(3.7) of the RW equation for any angular momentum. Special points need not lie on the imaginary axis, as is clear 
by considering V{x) ~ e~^^ cos(/ia;); however, only if Rew' = can w'^ be the eigenvalue of the SUSY generator, 
so this case needs to be studied most carefully. An example in which the special point index n > 1 is provided by 
V{x) ~ xe^'^^ with n — 2, with obvious generalizations. 

Physically, the divergence in gi at a special point is irrelevant: it merely comes from the normalization imposed on 
g. Thus, this divergence cancels against the Wronskian in the Green's function, the proper outgoing wave near uj = lo' 
in fact being x('^)- Since x(w) ^ [u — Lo')^e^'^^, the large part of x vanishes at the special point. Because x(^) is a 
solution of the KGE with eigenvalue o;^, this can only happen if the outgoing wave x('^')i ^^-id to be special outgoing, 
at the same time is incoming as well: x(<^') g{~^')- While this might be counterintuitive, it merely means that at 
these frequencies the exponential tail scatters so strongly that the true outgoing wave is completely different from the 
one in free space. All this can be verified beyond the Born approximation for the exactly solvable cases V{x) ~ e~^^ 
^ and the PT potential of Section 0. 

On the negative imaginary axis, the analytic continuation in general yields two different branches, g\{uj) = 
limexo .g(t^ - e) and g^{uj) = lim^io + e) = 9\{^)* ■ 

The type of discontinuity this induces is indicated for each lo' ~ ~i\^'\ by a second index m so that Ag{Lo) = 
gr{Lo) — g\{o-!) ^ {oj ~ uj')"^. The case in which there is no cut at all can be thought of as m = oo. Note that because of 
their normalization, the "large" parts of gi/^ cancel in Ag, which, also being a solution of the KGE, must thus be small 
if it does not vanish: Ag(Lo) oc g{—Lo), with a purely imaginary (frequency-dependent) constant of proportionality. At 
a special point, xi and Xr must be proportional, their prefactors having opposite phases, and the cut index is then 
defined by Ax{lo) ^ {lo ~ w')". 
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The above treatment might seem to put undue emphasis on exceptions. However, these cases will actually arise 
by performing SUSY transforms on a generic potential, as will become clear in the following. The statements in this 
subsection have obvious counterparts for /. 



3. SUSY transformation of one outgoing wave 



We first consider how the outgoing wave at one end is transformed under SUSY. The crucial feature to be derived 
is: even in the presence of tails, SUSY carries outgoing waves into outgoing waves unless it annihilates them. The 
latter exception can only occur if the eigenvalues of the outgoing wave and fi^ of the SUSY generator coincide. 



for Ag = iff (? oc $. The proof is in fact simple. By (1.7), SUSY carries solutions of the wave equation into 
solutions for the partner system. For Imw > 0, the operator A of ( |1.6| ) carries exponentially decreasing waves into 
exponentially decreasing waves since W is bounded. Since the frequency independent A maps analytic functions to 
analytic functions, clearly it commutes with the analytic continuation to the rest of the w-plane, completing the proof. 

Normalization of the transformed wave is straightforward, with the result that the outgoing wave in the partner 
system reads 

for generic w. In particular, this means that special points, branch cuts etc. are all conserved under SUSY if uj'^ ^ fi^. 
For = n^, one has to study the limiting behaviour of (Bl) more carefully. Writing ft — zLiK {K > 0), one has the 
following cases. 



(al) — K and —iK is not a special point. In this case $ oc g{iK), so that the zero in the denominator of (Bl) 
cancels, as it has to because g{uj) is necessarily r egula r for Imoj > 0. Hence, g{iK) = limj^^jif ^(w) = —iAgi{iK), 
with gk{iK) = kl~^dl^g{uj)\iK, cf. the notation of ( ^.13| ). The expression for giiK) is nonzero because giliK) g{iK), 
as follows from {H + K^)gi{iK) = 2iKg{iK) ^ 0. 

Let the original system have a discontinuity index m at frequency —iK. Frequency differentiation of the KGE 
yields {H + K'^)gk{—iK) = —2iK gk-i{—iK) + gk-2{—iK) on both branches, so that {H + K^)Agk{—iK) vanishes 
for k ~ m, but is nonzero for k = m + I. Thus, Agmi^iK) cx g{iK) (since Agm{—iK) is always exponentially 
small, cf. at the end of Section |B 2| ) and hence Agk{—iK) = for < fc < m (for k < m this follows simply from 
Agk{—iK) — 0). However, Agm+i{—iK) ^ so that m = m + 1. In particular, m > 1 so that aX lu — —iK, the cut 
discontinuity vanishes in the partner. 

It is trivial that $ is exponentially large, but is it actually outgoing? The question is answered by a calculation of 
the Wronskian 

j(g-\iK),Ag{-iK))^Q , (B2) 

for which one has to substitute ( |1.6D for A with W = —g'{iK)/g{iK) and use the KGE for g" . Thus, up to 
proportionality $ equals the unique outgoing wave g{—iK) = Ag\/-c{—iK)/2K. 

(a2) W+ = K and —iK is a special point with index n ~ \. The situation at = iK is identical to case (al). Near 
Lo — —iK, g{Lo) — x(tt')/(a' + iK) with %(— iif) c>c g{iK) (cf. Section |B 2| ), so that Ax{—iK) — and the special point 
is lifted in the partner system, with g{—iK) — Axi{—iK)/2K ^ 0, where the r.h.s. in general will be different on the 
two branches. In more detail, a calculation as for case (al) shows that now fh — m. 



The counterpart to (B2) is evaluated to be 



J(^g-H^K),Axli-^K)) = -2iKx{-^K)/gi^K) 

7^0 , (B3) 

so that in this case $ is mixed in the partner system fi^ . 

(a3) — K and —iK is a special point with index n > 2. The situation at lu — iK is identical to case (al). At 
to = —iK one calculates n = n — 1 with x(— iX) — Axi{—iK)/2K and fh = m. The function $ is exponentially 
increasing, hence certainly not (special) outgoing. All these properties are seen to be analogous to case (a2), except 
that now the frequency —iK remains special in the partner system. 
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(bl) VKf = —K, —iK is not a special point and $ cx g\{—iK) = gr{—iK) is outgoing. Near uj = iK the transfor- 
mation is completely regular and g{iK) = ~Ag{iK) /2K cx since $ is an exponentially decreasing solution with 
eigenvalue —K'^. 

Near uj = —iK, g{uj) = Ag{uj)/{iuj — K) ~iAgi{—iK) — g{~iK). This generalizes the situation of (S.17), where 
-zif is taken to be a double QNM. Setting A5(ti;) ~ (w+ii^T)™ (m > 1) leads to Ag{u;) = ~i{u;+iK)"^-^AAg,ni~iK)+ 
0[{uj + iK)"^] where AA.gm{—iK) is readily checked to be nonzero, so that m — m — \. 

(b2) W+ = —K, —iK is not a special point, and $ is mixed. Cf. case (bl): an exponentially increasing $ will 
always be mixed if there is a discontinuity across the cut, for in that case gi/i.{—iK) are not proportional to a real 
function. The situation at to = iK is identical to case (bl). 

Near lo = —iK, the zero in the denominator of (Bl) is not compensated, so —iK becomes a special point with index 
n = 1 in the partner system and x(— zfT) = —iAg{—iK). Note how the large parts of <& and g{—iK) coincide, leading 
to the annihilation of this part in Ag{—iK) and leaving only an exponentially decreasing function. One verifies that 
ifi = m for the discontinuity index. 

Note that, if V for instance has finite support, V has a tail ^ e~^^^ as in (2/7), while —iK is its o nly special 
point. Thus, at iuj/K = 2,3, .. . one has a cancellation of the kind called "miraculous" in Appendix B2, but which 
apparently is not exceptional in the context of SUSY. 

(b3) W+ = —K and —iK is a special point. The situation at = iK is identical to case (bl). Near lo = —iK, the 
pole in (Bl) leads to fi = n + 1 for the special point index with x{—iK) = —iAx{—iK), while rh = m. 

One sees that, for each i — 1, 2, 3, the partner of a system in the (az) category falls into case (hi) and vice versa. 
Thus, for each i, (a?) and (hi) transforms are each other's inverse, and inspection of the above cases confirms that 
their effects are exactly opposite. 

Again, the situation for x — » —oo is wholly analogous. 



4. Two outgoing waves 

We now combine the left and right boundary conditions to study modes. For this purpose, we introduce the two 
Wronskians 

J^{u;) = r{u;)g{u;)-f{u;)g'{u;) , (B4) 
MLu) = f'{-iu)g{iu)-fi-u;)g'iu) , (B5) 

which keep track of (quasi)normal modes and total-transmission modes respectively. That is, for ordinary points of / 
and g an Afth-order zero in Jq{uj—ujo) signifies an Afth-order QNM if ImajQ < 0, and an Afth-order NM if Imcjo > 0. 
In the latter case, this of course means that the only possible nonzero value is M = 1, in which case Rewo — 0. 
Similarly, an A/th-order zero in Jt(ci;=a;o) signifies an A/th-order TTMl at luq, or equivalently an A/th-order TTMr 
at — wq. 



While QNMs can coexist with NMs, or TTMlS with TTMrs, at the same eigenvalue uj^ (cf. Section |VII]| ), (Q)NMs 
and TTMs clearly exclude each other unless there are special points, at which the various notions coincide. 

Of course, if ujq is a special point of index n for / and n' for g, then (Q)NMs are given by zeros of (uj — ljq)"^" Jqi'-^) 
and TTMlS by those of (lu - wo)"'(w + wq)" Jt(t^)- 

For example, let g{u>) have a special point with index 1 at w = loq. Then at the eigenvalue lOq one can have: («) no 
mode, or (ii) an NM and a TTMl of order 1 + A: (fc = 0, 1, . . .) or (Hi) a QNM of order 1 + i and a TTMr of order 
1 -|- A: (A:, ^ = 0, 1, . . .), with at most one of k, I nonzero. As a second example, if wq is an n = 2 special point for /, 
then there can at most be a first-order TTMr at that frequency, since now the notions of TTMr and NM coincide 
up to second order and higher-order NMs cannot occur. 

If at least one of the functions / and g has a cut on the negative imaginary axis, then Jq will have such a cut as well 
(similar to the cut in the Green's function). On the other hand, clearly if there is a branch cut in / (g) then Jt will 
have a cut on the positive (negative) imaginary axis. It follows from the symmetries satisfied by / and g individually, 
that if is a zero- mode at frequency loq on the left branch of Jq, then 0^. — (f>^ is a zero- mode (of the same order 
M) on the right branch. As a consequence, if fi{ujQ) = /r(wo) (for instance if / has no branch cut at all) then 5(^0) 
should be proportional to a real function; hence, if loq is an ordinary point for g then Ag vanishes there (mR > 1). 
For TTMs the situation is simpler, since only one cut at a time is involved: if, say, there is a TTMl at ujq, then 
max(mR, wr) > 1. Generalizations to higher-order modes can be derived, but are not given here. 
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5. SUSY transformation of two outgoing waves 



With 



the 



SIX 



possibihties 



of 



Section 



BS 



for 



both 



X —oo and x ^ oo, and with the rich variety of possible mode structures in the original system outlined in 
Section B 4 , a classification of the possible SUSY transformations may seem daunting even before considering the 
subtleties of = 0. However, most of the information is contained in the transformation properties of the Wronskians 
(B4) and (B5). These are readily derived by substituting (Bl) and its counterpart for / into (84) for Jq and ( p35| ) for 
Jt, using (1^) for A and the KGE for /" and g" . One obtains 



{uj + iW-){LJ ~iW+) 

[uo + iK){uj - iK) 
{-LJ + iW-){uj -iW+) 



(B6) 
(B7) 



Combining these with the transformation properties of the individual functions / and g derived in Section B3, one 
finds that also the mode structure is entirely conserved if 7^ leaving only the case — il^ to be studied. 

For instance, if W± = ±K (an NM), then (|B^ ) shows that in Jq a zero at iK gets cancelled. This is a much more 
meaningful statement than just noting that = for the NM SUSY generator <&, which in itself says nothing about 
the situation in the partner system. At the same time, an extra zero at —iK appears. However, this latter zero must 
be interpreted carefully. If —iK is an ordinary point for both / and g, it will remain so in the partner and the zero 
signifies that a QNM at —iK is created (or its order increased). If —iK is special for either / or g, then the extra 



zero merely reduces n by one according to case (a2) or (a3) in Section B3, and the order (possibly zero) of the QNM 
at —iK is not altered, consistent with $ not being outgoing on the 'special' side. Finally, if both / and g are special 
at —iK then the decrease of the associated indices by itself would reduce the pole in Jq by two orders, so that in this 
case (B6) shows that the NM-generated SUSY transform reduces the QNM index at —iK by one. In fact this should 
not be surprising, for if —iK is two-sided special then at least to first order the QNM agrees with the NM, which gets 
cancelled in the partner; an example is the inverse of the transform considered in Section [VE| . 

Continuing with the example, one sees that Jt merely changes sign. Again, depending on the properties of / and 
g, this can mean no change in the mode situation (if —iK is ordinary), destruction of a TTMl (if —iK is special for 
g), destruction of a TTMr (if —iK is special for /), or both of the latter. 

Generalizing from this example, one sees that the material of Section [B 3| together with ( p6| ) and (B7) leads to the 
following statement: if —iK is an ordinary point for both the original and the partner system (which in particular 
implies that $ is purely incoming or outgoing on both sides), the only effect of SUSY on the mode structure consists 
of the destruction of the mode 3> and the creation of the mode in both cases possibly by changing their order by 
one. Hence, in this case the situation remains as in Table 1 even in the presence of arbitrary potential tails. 

It may seem that (B6) holds the possibility of creating a second-order NM by starting from a system with an NM at 
iK, and performing a SUSY transform using a generator for which W± = ^K. However, this can never occur because 
at the eigenvalue —K^, any real solution except possibly the NM itself must have a node and hence cannot be used as a 
SUSY generator. For a proof which is also valid at K — 0, assume that </) is the NM with 0(a;^oo) > for definiteness, 
while $ 9^ (/) is a positive-definite solution at the same eigenvalue. As x 00, $/0 ^ 00 so ($/</>)' = J(^, (t')/4>^ > 0- 
Since the Wronskian is position independent, (<&/</>)' > also for x —00, which is only possible if </) < there, 
implying that the function has at least one node. At the rightmost node, J($, 0) = -cj)'^ < 0, and the contradiction 
with the position independence of J means that the bmction $ of the assumption does not exist. 



6. SUSY at f7 = 



Throughout this appendix and the main text, it has been clear that there is an interplay between the outgoing waves 
and possible modes at iK and —iK. Hence, special attention should be paid to the case when these two frequencies 
coincide, i.e., f2 = 0. For the KGE at a; = 0, it is still true that for x —^ 00, there is one unique (up to a prefactor) 
small solution. This same solution is also the outgoing wave x(^=0), as follows by analytically continuing the small 
solution g{uj) from the upper half w-plane to the origin (the same argument can not be used to conclude that x('^=0) 
should be asymptotically large by performing the continuation from the lower half plane, since at any frequency large 
functions are not unique, and the limiting member of a sequence of them could well be asymptotically small). The 
issue of special points does not arise at zero frequency, and neither does that of discontinuity indices. 
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For = 0, the six cases of Section B3 reduce to two: for x oo, $ can be either (a) outgoing, i.e. small, or 
(b) large. However, the properties of the transformation depend sensitively on the potential tail. Here we study 
V{x) k/x", with a> 2 for the practically important case a = 2 we take k — + 1) > —1/4, with v > —1/2 
without loss of generality. Tails which decay faster than algebraically (including no tail at all) are easy to handle, and 
can be thought of as a = oo. 

For a > 2, all terms in the Born approximation exist down to cj = 0. To first order p^ , g{0,x) — 
l + ra^""/ (a— l)(a— 2), so in case (a) one has W{x) — Kx^^"/(a— 1), i.e., SUSY reverses the sign of the potential tail. 
In the partner system, ^(0) = hm^^^o Ag{uj)/iu; = —idujAg{uj)\f^ oc $; only if a > 3 can one write ^(0) = —iAgi{0) oc l> 
for this, since g{oj) itself is not differentiable at = for 2 < a < 3. The proportionality Agi{—iK) cx $ has not been 



encountered before in Section B3; Eq. ( p33| ) shows why. In case (b), ^ ^ x + 0{x^~") + 0{x^), so W{x) ^ —l/x, 
V{x) ~ 2/x^, and g{uj) — x['-^)/^ with $ (x x(0) — —iAg[Q) ~ i/x. Thus, for these potentials at O = 0, any choice of 
$ leads to a small $, in contrast to the situation at finite SI; of course, the partner potential depends on this choice. 
If a = 2, one has the exact solution giuj, x) ~ e^'^^^^^'^ 1'^ ^ -klox l2B.'^^]^^i^(tjJx^ ||l^ with a Hankel function, i.e., 

g{uj,x) ^ [T{l>+^)/^y¥\{2i/uJxY (~ {i— I) In ujy/ujxjn for u = — |) near uj = 0. In case (a), the generator thus is 
$ — and the transformation can be summarized by = i/ — 1 (— t^) for > ^ (— | < ^ In case (b), one has 
<i> ~ x"^^ and z> = + 1 |48|. Thus, one sees that for — i < < i, always $ oc x(0) (indeed, a large $ would give 

i-^ i> + 1 in the inverse transform, which for any i> is incompatible with < i) as for a > 2. Such potentials, which 
at = behave essentially the same as one with finite support, have weak tails. On the other hand, for strong tails 
(a = 2, > i) $ is large (small) if $ is small (large), as is the case for J7 ^ 0. 

In line with usage for other frequencies, a; = is said to be a mode frequency if x{'^=0) (/)(w=0), where </> is the 
counterpart of x for waves outgoing to the left. Like any NM, such a mode is necessarily of first order, and there is 
no distinction between different types of modes. Apart from possible overall prefactors 4>{uj)/ f{uj) and x(a;)/(7(a;) due 
to the tails, a mode is associated with a second-order zero in Jq/t because the eige nvalue in the KGE is w^. 

If the original system has a mode at w = 0, the proof at the end of Section B5 implies that (if it is nodeless) the 
mode itself is the only candidate for a SUSY generator, i.e., case (a) applies on both sides. If V has two weak tails, $ 
will again be a mode; combination with the above statements on the transformation of one potential tail shows that 
the class of weak-tail potentials with a mode at w = is clo sed under SUSY (since transforms at 7^ obviously 
cannot generate strong tails), as exemplified in Section VE by some of the potentials SUSY-equivalent to the free 
field. If V has a weak tail on one side, $ is outgoing on that side only so that there is no mode at w = 0. liV has 
two strong tails, $ is large on both sides. Moreover, according to the above in that case one has t'L/R = ^'l/r ~ 1 for 



the tail amplitudes, so that f {uj)g{uj) / f (w)g(ijj) ~ w^; since ( |B6| ) and (B7) show that up to a sign Jq/t are invariant 
under SUSY, the compensating factor cj~^ must come from the disappearance of the mode. 

If the original system has no mode at w = then the situation is simpler, even though several types of SUSY may 
be possible. If <E> is outgoing on both sides (which in particular will always be the case if V has two weak tails) then 
it is obviously a mode itself. It is not iff <& is small on only one side and if V has a strong tail on that side. Then 
is small on only the other side, i.e., if $ is not a mode then V has no mode at all at w = 0. 

To conclude, let us mention that modes, special points etc. at 7^ = are not affected by SUSY, as is the case 
for other generator eigenvalues. 



APPENDIX C: NODES IN QNM WAVEFUNCTIONS 



Nodeless eigenstates play a special role in SUSY: they are candidates for the SUSY generator. It is therefore 
useful to highlight the differences between NMs and QNMs in this regard, especially to contrast with the well-known 
property that there can be only one nodeless NM. 

First of all, we show that for QNMs not lying on the imaginary axis (i.e.. Re a; 7^ 0), there can be at most one 
node or antinode. This is not surprising: since the eigenvalue is complex, the wavefunction has a changing phase, 
and it would be "unlikely" that the real part and the imaginary part (or their derivatives) would vanish together. To 
prove this, we take the Schrodinger point of view, so that the eigenvalue is A = with a nonzero imaginary part. 
Now consider a time-dependent QNM and suppose that it has nodes or antinodes at two points xi,X2- At these two 
points, the current 

J = z[0*(9,0)-(5,0*)(/.] (CI) 
vanishes. Then, applying the conservation law in integral form to the interval [xi, ^2], we find that the total probablity 
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in this interval is constant in time. Yet the wavefunction is either growing or decaying, since Im A 7^ 0, which is 
therefore a contradiction. 

From the perspective of SUSY it is unfortunate that the above proof excludes the crucial imaginary axis. However, 
on that axis the statement remains valid for repulsive potentials, or more generally for potentials which are so weakly 
attractive that V — oj'^ is positive definite. Namely, let be a solution with two (anti)nodes. By taking the real or 
imaginary part if necessary, (jj can be arranged to be real. Now between two nodes (j) would have an extremum, i.e. 

< which is incompatible with the KGE. Similarly, an antinode can only be a global maximum or minimum, 
precluding the presence of any other nodes or antinodes. Note that for repulsive V, the KGE can be mapped to the 
wave equation (Appendix the mapping conserves nodes but not necessarily antinodes) , for which a proof using the 
current J (cf. (pl|)) does not exclude the imaginary axis. 

Thus, QNMs (at least those off the imaginary axis for attractive V) can have at most one node or antinode. For 
symmetric potentials, in the even sector x = is already an antinode, so there can be no nodes anywhere. 

For zero modes, i.e., QNMs with Rew = 0, nodes are more "likely": the eigenvalue is real and (if there is no branch 
cut) the wavefunction has a constant phase, say purely real; it therefore requires only one condition, rather than 
two, to have a node. Nevertheless, in contrast to the conservative case, the proof that there is at most one nodeless 
eigenstate can be bypassed. 

The interlacing nodal structure of NM eigenfunctions follows from well-established Sturm-Liouville theory. For the 
present purpose, we do not need the full apparatus. Consider, for simplicity, a finite interval [~a, a] and suppose there 
are two distinct nodeless eigenfunctions 0i , 02 , both of which can be chosen to be real. Then they can both be chosen 
to be positive, which violates the orthogonality condition for NMs 



Mx)Mx)dx = . (C2) 
We can attempt to transplant the argument to QNMs. For zero modes, again the wavefunctions have constant 



phases and can be chosen to be real, and if they are nodeless, positive definite. However, the analog of ( |C2| ) for two 
eigenfunctions with eigenvalues uJi — ~i"fi is 



-(71+72) / (t>i{x)(l)2ix) dx 

J —a 

+ [M~a)M~a) + Ma)Ma)] = ■ (C3) 

Note in particular the signs of the two terms. With 7^ > 0, this condition does not preclude both eigenfunctions from 
being positive definite. 

Thus, we have three different situations, (a) For NMs, there can be only one nodeless state, (b) For QNMs with 
Rew = 0, there could be more than one state with no node, (c) For QNMs with Rew ^ 0, each eigenfunction can 
have at most one node or antinode, and for symmetric potentials, every even eigenfunction is nodeless. 



Case (b) in particular opens up the possibility of multiple SUSY transformations, as exemplified in Sections HIB 
and 



APPENDIX D: SELF-REPLICATING POTENTIALS 



In Section |V C| we saw that the PT potential is self-replicating under SUSY, in the sense that the partner is the 
same potential with a different strength: 



V{x) = aV{x) 



The condition (Dl) is a special case of shape-invariant potentials described by a family U{a,x) [[49|, with 



V{x) = U{a,x) 

V{x) = U{a,x) + c{a) 



(Dl) 



(D2) 



Here we focus on the more restrictive sense (Dl) and show that the PT potential is uniquely determined by this 
condition. From (Dl) and (2.1), and with Q = —iK being purely imaginary, we find 



W'{x) = (3 [W{xf --K^] , P = 



a— 1 
a-fl 



(D3) 



26 



All nonsingular nonconstant solutions of (|D3|) differ merely by x-translation. The unique antisymmetric one reads 

(D4) 



W{x) = -KtSinhK/3x , 



Putting this into (2J), we then find 



V{x) = if3-l)K^ sech^ K(3x 
V(x) = i-P-l)K^sech^ K^x 



(D5) 



That is, V and V are of PT form. 
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potential so that its small solution x('^=0) is $ = instead of a;~^/^. 
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29 



Type 


$ 


$ 


NM 


QNM 


TTMl 


TTMr 


1 


DD 


II 


-1 


+ 1 








2 


II 


DD 


+1 


-1 








3a 


DI 


ID 








-1 


+1 


3b 


ID 


DI 








+1 


-1 


4 


more complicated, see Appendix ^ 



Table 1: The boundary conditions for the deleted state <& and the new state $, and the changes in the number of 
NMs, QNMs, TTMlS and TTMrs in the SUSY transformation. 
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0.0 0.5 1.0 

Fig. 1(a); A squarc-wcU potential V (solid line) defined by (p.ll) with Vq 



-20 and a — 1, and its SUSY partner potential V (broken line). 



A 



-10 -5 5 10 

Fig. 1(b): The complex u;-plane showing the NMs and QNMs common to both potentials (crosses); the mode present only in V (circle), which 
corresponds to the generator <1>; and the mode present only in V (triangle), which corresponds to ^ — A^. 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Fig. 2(a): A square-barrier potential V (solid line) defined by ( |j.l| ) with Vb — 0.16 and a — 1, and its SUSY partner potential V (broken line). 
The SUSY transformation is construeted by using the state at uJi — — O.lSli (cirele in Fig. 2(b)) as the generator. 
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Fig. 2(b): The complex a;-plane showing the QNMs common to both potentials (crosses); the mode present only in V (cirele), which corresponds 
to the generator and the mode present only in V (triangle), which corresponds to ^ — A^. 
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"O.O 0.2 0.4 0.6 0.8 1.0 1.2 

Fig. 3(a): A squarc-bariier potential V (solid line) defined by ( |3. l| ) with Vq = 0.16 and a — 1, and its SUSY partner potential V (broken line). The 
SUSY transformation is constructed by using the state at UJ2 — —2.5002 (circle in Fig. 3(b)) as the generator. Inset shows one portion enlarged. 
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Fig. 3(b): The complex a;-planc showing the QNMs common to both potentials (crosses); the mode present only in V (circle), which corresponds 
to the generator ^; and the mode present only in V (triangle), which corresponds to ^> — A^. 
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-1.2 -0.8 -0.4 0.0 0.4 0.8 1.2 

Fig. 4(a): A multi-step potential V (solid line) defined by ( [j.ol ) with Vq — 1.0, Vi — —10.0, a — 1.0 and b — 0.1, and its SUSY partner potential 
V (broken line). 



-10 -5 5 10 

Fig. 4(b): The eomplex oj-plane showing the NM and QNMs common to both potentials (erosscs). The square indieates a TTMl and a TTMr in 
V, and a doubled TTMr in V. 
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Fig. 6 (a): A square-barrier potential V {solid line) defined by ( [j. l| ) with Vo — 0.16 and a — 1, and its SUSY partner potential V (broken line) 

generated by a Type 4 transformation, using K — 3, c — 1.0 and d — —0.829. Inset shows one portion enlarged. 



-6 6 

Fig. 6(b): The complex [A;-plane showing the QNMs common to both potentials (crosses); the mixed mode $ with boundary condition MM in V, 
used as the generator (square). 
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Fig. 7(a): The truncated PT potential (5.11) with V — 3/16, 6—1, and a — 2, and its SUSY partner potential V (broken line) generated by , 
Type 2 transformation, using the state at £7 — —0.224*. 
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Fig. 7(b): The complex a;-plane showing the QNMs common to both potentials (crosses); the mode present only in V (circle), which corresponds 
to the generator ^; and the mode present only in V (triangle), which corresponds to ^> — A^. 
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Fig. 9: The distribution of QNMs of a Schwarzschild black hole for I = s = 2. The mode on the imaginary axis 2mw = —4i is not a QNM, but the 
special mode. 
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